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The close-limit method has given approximations in excellent agreement with those of numerical 
relativity for collisions of equal mass black holes. We consider here colliding holes with unequal 
mass, for which numerical relativity results are not available. We try to ask two questions: (i) 
Can we get approximate answers to astrophysical questions (ideal mass ratio for energy production, 
maximum recoil velocity, etc.), and (ii) can we better understand the limitations of approximation 
methods. There is some success in answering the first type of question, but more with the second, 
especially in connection with the issue of measures of the intrinsic mass of the colliding holes, and 
of the range of validity of the method. 
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I. INTRODUCTION AND OVERVIEW 

> . 

Recently much attention has been given to the problem of the collision of two black holes. There are two main 
^ \ \ reasons for such interest. First, black hole collisions provide one of the most powerful and interesting sources for 
possible detection by gravitational wave observatories W. Second, black hole collisions are being intensively studied 
numerically by using supercomputers to evolve Einstein's equations from initial data representing two holes [Q . 
t-H . The enormous difficulty of solving the equations numerically provides motivation for approximation methods that, 
with little effort, can give guidance to what results can be expected, what cases are interesting for full numerical 
study, etc. One method that has proved surprisingly successful in a range of tests is the close-limit approximation 
(CLAP). This method applies to initial value data representing holes which are initially close to each other. If the 
holes are close enough then the horizon will initially surround both holes, and the spacetime outside the horizon can be 
I ■ considered to be a single perturbed hole. Strong nonsphericity inside the initial horizon cannot affect the evolution of 
\Q j the fields outside, and hence does not influence the generation of outgoing radiation. The nonspherical perturbations 
On ■ outside the horizon can be analyzed with the well developed techniques of linear perturbation theory 

' The method, applied to numerically or analytically generated initial data was discussed by Abrahams and Price 
and has been applied to: (i) simple analytic initial data for the head-on collision of momentarily stationary black 
holes Q-^], (ii) numerically generated initial value data for holes which are initially moving towards each other [Q, 
j^Tjl (iii) analytic initial value solutions for holes which are initially moving slowly towards each other (iv) analytic 
initial data for holes which have opposite inital angular momentum and are initially momentarily stationary ||. 
Where comparisons with full numerical results are available (all of the applications above, except the last) the results 
of the CLAP method are found to be remarkably successful, even when initial conditions would seem to violate the 
assumptions underlying the approximations. This success holds the promise of giving easy approximate answers about 
9^ , black hole collisions. 

All applications listed above have been to collisions of holes with equal masses. (Those are essentially the only cases 
for which comparisons have been available with fully numerical computations.) Here we consider collisions of unequal 
mass holes, though comparisons with numerical results cannot yet be made. There are two reasons we do this. First, 
there are some interesting questions to which it is better to have a very rough answer, even an uncertain answer 
(if the nature of the uncertainty is kept in mind), than no answer at all. It is interesting to see whether there are 
any surprises in predicted radiation efficiency, in dependence on details of initial data, etc. One interesting question 
applies only when the masses are unequal: can gravitational radiation from the collision contain a significant amount 
of linear momentum, so that the hole that forms will recoil with an astrophysically significant proper velocity? 

The second motivation is to look at the range of validity of CLAP calculations. The ultimate index of validity (short 
of full numerical relativity comparisons) is second-order perturbation calculations, and work on this technique is well 
underway Jio| , pd| . These necessary calculations, however, are quite difficult (though far easier than full numerical 
computations) and it is useful to look at whether simple guidelines exist. Our approach here is simultaneously to use 
the CLAP method to look for astrophysical answers, and to use the examples to gain deeper understanding of the 
method. 

The nature of the questions being asked justifies avoiding unnecessary complications and using the simplest initial 
data sets applicable. We therefore limit our attention here to nonspinning holes which start from rest and undergo a 
head-on collision. More specifically, we limit our investigation to two momentarily stationary initial solutions. Both 
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initial geometries are conformally flat, and therefore are completely determined by a conformal factor <f> 4 , where $ 
solves the flat spacetime Laplace equation jl2| . The simplest solution is that of Brill and Lindquist (hereafter 
"BL") in which <f> has the form of the gravitational potential for two Newtonian point masses. Another solution, that 
of Misner and Lindquist [^4|-^6| (hereafter "ML" ) is more complicated but has the very useful feature of an easily 
located minimal area of the Einstein-Rosen bridges in the initial geometry. It is also the initial geometry that has 
been used in almost all numerical relativity studies of collisions of black holes. For these two sets of initial data we 
look at radiated energy as a function of how close the initial holes are. To characterize the separation of the holes, we 
use in all cases the proper distance L, along the symmetry axis, between the apparent horizons of each of the holes. 

An important issue that arises is how to characterize the "bare mass" of each of the throats, i.e., how to assign an 
intrinsic mass to each of the holes participating in the collision, a mass unaffected (in some sense) by the presence of 
other nearby sources of gravitation. There is a fairly natural choice of bare mass for the BL initial data, but not for 
the ML initial data. We consider three candidates as bare mass measures in the ML case. One of the more interesting 
conclusions of this work is the importance of the choice of bare mass, and the unphysical consequences of the wrong 
choice. 

In Sec. II we start by describing the BL and ML solutions and the choices that can be made for bare mass. We 
then describe the CLAP method. CLAP estimates of radiated energy are presented in Sec. Ill where it is seen that 
the choice of bare mass of one ML hole governs even the qualitative nature of some answers. In Sec. IV results are 
given for the recoil velocity of the final hole formed, due to the emission of radiation. In the results of Sees. Ill and 
IV, the usefulness of the results depend heavily on the range of validity of the CLAP. In Sec. V we look at a simple 
criterion for when linearized theory should be applicable. We find that the validity of this criterion, which seems 
useful in the case of equal mass holes, does not seem to extend to collisions in which the ratio of hole masses is very 
small. Conclusions are presented in Sec. VI. Details of several calculations are presented in three appendices. 



II. APPROXIMATION METHOD AND INITIAL GEOMETRIES 



Initial value solutions 



For simplicity we limit attention to time symmetric initial data. There is then an initial hypersurface on which the 
extrinsic curvature is zero so that initial value data for Einstein's equations consists only of the three geometry of 
this hypersurface. Also for simplicity, we limit our consideration to three geometries that are conformally flat. This 
simplification is neither selected for or against by any strong physical argument, but it leads to a very convenient 
mathematical description, and is therefore used, e.g., in most numerical relativity work. 

We write the conformally flat axisymmetric three metric in the form 

ds 2 = $ 4 (i?, 9){dR 2 + , (I) 

where dQ 2 = d8 2 + dip 2 sin 2 9, and R, 9, <p are spherical coordinates in the flat conformal space. The Einstein initial 
value equations then turn out to require only that $ obey Laplace's equation 

V 2 $ = , (2) 

where V 2 is the Laplacian with respect to the flat metric dR 2 + R 2 dH, 2 . 

It is necessary to find solutions corresponding to two initially static black holes. To investigate the possible sensitivity 
of radiation to the details of the initial value solution we consider two different solutions of (||) . Brill and Lindquist 
investigated the solution to (^) that is most immediately apparent Jl3| , the solution with the form of the Newtonian 
gravitational potential of two mass points: 

$ = 1 + ^L + ^^. (3) 
\R-Ri\ \R-R 2 \ 

Here all the vectors and their norms are defined in the flat three dimensional space, with R being the position of an 
arbitrary point in such space and Ri is the position of a point in the flat space representing hole i. Though (^) suggests 
two pointlike solutions, there are in fact no physical singularities corresponding to these points. Rather the three 
geometry near each "point" can be extended, through a throat, out to an asymptotically flat space (an Einstein-Rosen 
bridge). The complete BL three geometry, then, includes three asymptotically flat regions. One is a region with two 
throats connected, the region of "our universe." The other two regions contain one throat each corresponding to mass 
1, and to mass 2. (These distinctions are meaningful only when the two throats are well separated in "our" universe; 
see p|). 
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The BL topology has a useful practical feature. Since hole 1 (for example) has its own asymptotically flat region we 
can infer a "bare mass" mi for hole 1 from the metric at large distances from the throat. In this way we can compute 
bare masses rrii for each of the individual holes. (In the limit that the holes are very far apart, the bare masses are 
related to the parameters on in (||) by on w TOj/2, where we use c = G = 1 units.) We can also compute the total 
mass M in "our" universe, the ADM mass of the two-hole spacetime. 

Another solution to (|J), corresponding to two throats, of the form 



can be constructed by placing the "masses" a n and b n at the locations d n and e n in the flat conformal three geometry 
in a manner similar to the placement of electrical charges in the problem of finding the electric potential distribution 
in the region outside two charged conducting spheres. With an infinite set of image "masses" it is possible to construct 
a three geometry with two throats which open into two identical asymptotically flat regions. The isometry between 
the two regions takes the form of a reflection through spheres at the minimal neck of each of the throats. The details 
of the symmetrization procedure are given in Appendix B (see also Lindquist jig] ). 

For our purposes, the disadvantage of ML solutions is that there is no immediate meaning that can be given to 
"bare mass." We consider three candidates for bare mass in ML solutions. One candidate is the bare mass suggested 
by Lindquist in his study of symmetrized initial value solutions [fj6| . The "Lindquist mass" sums the bare masses 
of all the images associated with one of the throats (a divergent sum) then subtracts a Newtonian expression for 
the binding energy due to the interaction of those images with each other (another divergent sum). This definition 
of bare mass seems not to have been used in recent work on the problem. Much more commonly cited is a rather 
straightforward measure of the bare mass that we shall call the "area-mass." One takes the area Ai of the minimal 
throat on the initial hypersurface and computes a bare mass from it as if it were an isolated Schwarzschild horizon: 
mi = Ai/16ir. A third candidate for the bare mass is the "Penrose mass" [T^-[l9| . This is a quasi- local definition of 
mass interior to a 2-sphere that starts with a method, valid in linearized theory, of extracting mass information from 
the Weyl tensor. This method is then formally converted to curved spacetime. The method, cannot be applied to all 
spacetimes, but it is always applicable to axisymmetric cases. Some details of the computation of the Penrose mass 
are given in Appendix C. 

FIG. 1. For two holes with fixed Penrose mass rriip, in the ratio m2p/rri\p = 3/7, the dependence of the Lindquist mass 
rriiL on separation is shown. 

Figure [l] illustrates the issue of choice of bare mass in the ML case with an example in which the ratio of Penrose 
bare masses of the two holes is 7:3. For each hole, the ratio of Lindquist mass to Penrose mass is shown as a function 
of separation. For both the heavy and the light holes, the Lindquist mass decreases relative to the Penrose mass, but 
the decrease is more dramatic for the heavy hole. If the statement "put the same two holes at different separation" 
means putting holes of fixed Penrose mass at different separations, then these "same" two holes at small separation 
have very different Lindquist masses, and a very different ratio of Lindquist masses. A similar comparison between 
Lindquist mass and area-mass tells a very different story. If the m 1 /m 2 ratio of area-mass is 7:3 then, to considerable 
accuracy, the ratio of Lindquist masses is also 7:3. There is not a strict equivalence of the two ratios. At very small 
separation there is a small deviation in the ratios. More important, the value of the Lindquist mass and that of the 
area-mass for a given hole are not the same when separations are small, but even here the effects are small. This 
is shown in Fig. g for two equal mass hole. At small distances (L/M less than around 0.5) the Lindquist mass is 
discernibly larger than the area-mass. By comparison, the difference in Penrose mass and area-mass is 20 times larger. 
In our energy and recoil estimates in the next two sections, no difference can be seen in the results depending on 
which mass ratio is held constant, the Lindquist mass or the area-mass. We shall, therefore, present only the latter. 

FIG. 2. Difference between Lindquist and area mass as a function of separation when the holes are equal. Here M is the 
total ADM mass of the spacetime. 

FIG. 3. For the bare mass U12 of a hole, the ratio of the Penrose measure of bare mass (m^p) to the Lindquist measure 
(wi2l) is shown as a function of the ratio of (Lindquist) bare masses for two different separations L/M (where M is the total 
ADM mass of the spacetime). The result shows that in the extreme ratio limit different measures of bare mass agree. 

In the limit of extreme mass ratio (as well as in the limit of large separation) the ambiguities of bare mass disappear. 
When the limit 1112/ mi of the bare masses becomes very small, the "particle" limit for m2, physical intuition suggests 




(4) 



3 



that any reasonable definition of bare mass must agree with the proper mass of a point particle perturbing the 
spacetime, and therefore all reasonable definitions of bare mass will yield the same result. This is illustrated in Fig. 
which shows the ratio of the Penrose and Lindquist measures of m-i , the bare mass of the less massive hole, as the 
ratio of bare masses 7712/TO1 goes to zero. It is clear that in this limit also all reasonable measures of mi, the bare 
mass of the more massive hole, will agree with each other and with the ADM mass of the spacetime. As a check we 
have computed the ratio of Penrose to Lindquist values of mi as mijm\ decreases and have verified that the ratio 
goes to unity (although not as quickly as in Fig. |J). 



The close approximation 

To apply the CLAP to the momentarily stationary axisymmetric space we must rewrite ([[]) in the form of a 
t constant slice of a Schwarzschild spacetime. This requires mapping the geometry of ([!]) onto a set of Schwarzschild- 
like coordinates r, 9, ip. We do this by transforming the radial coordinate as if R were the isotropic radial coordinate 
for the Schwarzschild geometry: 



R= (Vr + vV - 2M) 2 /4 , (5) 
where M is the total mass of the spacetime. This transformation puts the spatial metric in the form 

cfe 2 = W d f +r 2 <m 2 ) , (6) 



where 



1 - 2M/r 



(7) 

1 + M/2R ' v ' 



The function $ must satisfy (||) and will have singularities at the coordinate locations Ri of "mass points" (BL) or 
"images" (ML). Since $ must approach the Schwarzschild spatial metric $ — * 1 + 2M/R as R — > 00, we can expand 
for R > max (Ri), in Legendre polynomials Pi as 



2M 00 / M\ e+1 
$ = 1 + — +J> - P e (co S 0). (8) 

We can eliminate the £ = 1 term from (||) by appropriate choice of coordinate origin. The expression for J- ", therefore, 
can be put into the form 

For an initial value solution representing two black holes, the coefficients 7^ will contain a parameter describing 
the separation of the holes. As the separation goes to zero, the geometry approaches that of a t = constant slice 
of the Schwarzschild geometry, and the 7^ must therefore all approach zero. The close approximation consists of 
treating the separation of the holes as a perturbation parameter and in (|^) keeping only the terms lowest order in the 
separation. When we take the fourth power of J 7 the mixing of the Pe gives us, in principle, very complicated mixtures 
of contributions of different I for each final multipole of F 4 . The result in practice is much simpler. For both the BL 
and ML solutions (in fact for any conformally flat metric ([!]) for which (||) is used) the j£ coefficients increase in the 
order of the perturbation as £ increases and, as a result, for each £-pole of T the contribution of lowest perturbative 
order is due only to the 7^ term, and the conformal factor can be written 



M/2R 

' 1=1 

Equations 



^ 4 = 1+ i^m7^5>U Pe{cos9) - (10) 



([j]) , and ( [l0| ) describe an initial perturbation of the Schwarzschild geometry. Following the prescription 
of Moncrief ~pd| , we can use this perturbed initial solution to give us, for each £, the initial value of a "Zerilli function" 
y ipi. (We use the normalization of Abrahams and Price H, our ipg being their Q~\ .) 
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This Zerilli function satisfies a simple wave equation and from its initial value it is simple to find the time evolved 
function i/;(t,r). Once this is known we can compute the gravitational energy radiated by the system of two black 
holes from 



E 
2M 
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It can easily be shown that the Zerilli functions corresponding to two different solutions A and B of differ only in 
amplitude, i.e. they are related through 



This means that if we know ipi (r, t) for a given solution of (Q) , say the Misner equal mass solution, we can compute 
the energy (Sec. Ill) and the recoil velocity (Sec. IV) for any other conformally flat solution, directly. 



We give here the radiated energy corresponding to a range of head on collisions. To choose a particular black hole 
collision, we must first say whether the initial data is BL or ML. Second, we must specify a ratio of hole bare masses, 
and third we must specify L the separation of the holes, i.e., the proper distance between the apparent horizons. Our 
results are plotted with the radiated energy and the separation L normalized by the sum of the bare masses, so that 
the results can be interpreted as showing the change of radiation with separation when the "same" two holes are 
moved closer together. 

For ML holes we consider several different definitions of bare mass, and in Fig. || show the results for the Penrose 

FIG. 4. For two ML holes, with bare mass ratio F = milm\, the radiated energy is shown as a function of the proper 
separation L between the apparent horizons. Results are shown for both the Penrose bare mass definition and the Lindquist 
definition. For extreme mass ratios the choice of bare mass is irrelevant, bnt for F = 0.75 a qualitatively different dependence 
on separation shows up. 

definition and the Lindquist (cquivalently, the area) meaning of bare mass. To avoid cluttering the plot the results 
for F = 1 have been omitted. These F = 1 results show little difference depending on the choice of bare mass 
definition. The F = 1 curves are only slightly different for Penrose and Lindquist mass, and both are qualitatively 
similar to the results for F = 0.75, with Penrose mass. 

The results in Fig. [| for F = 0.01 show the expected independence of bare mass definition. But they fail to show 
another important feature that has appeared in a study of the particle limit (22). For particles that fall, starting from 
rest at distance L from a Schwarzschild hole, here is an anomalous decrease of radiation with decreasing L. This is 
found only in the range L = 5M to 7M, where M is the total ADM mass of the spacetime. This unexpected reversal 
of the general trend is a manifestation of the role played by the "curvature potential" (due to strong field effects near 
the hole) which peaks at around L = 3M. 

For our F = 0.01 results the ADM mass should be negligibly different from mi + m^- So the plot of radiated energy 
should show this effect in Fig. [|. The fact that it does not can be ascribed to failure of the close limit. We know in 
fact that this is the case from the study of the particle limit |2^] . For extreme mass ratios the close limit will work 
only for L less than about 2M, and when it applies the radiation decreases with decreasing L. 

For equal mass holes there is no anomalous region. Radiated energy is a monotonic function of initial separation. 
This can be interpreted to mean that when the smaller hole is not at the particle limit it is "too large" to be situated 
in the narrow range of radii for which the anomalous energy-separation behavior occurs. As we decrease the bare 
mass ratio F from unity, there must be some value at which the anomalous region appears. The results in Fig. ^are 
interesting in connection with this. For F = 1 it is well established that close limit results are reasonably accurate 
for values of L/M around 4. It would suggest that the anomalous behavior seen at rather small L for the F = 0.75 
is not an artifact of the close limit. The dramatic low-L bump in that curve, of course, has no counterpart in the 
corresponding Penrose curve, which suggests that what we are seeing is an effect due to an anomalous choice of bare 
mass. 

FIG. 5. For two BL holes, with bare mass ratio F = 1112/ mi, the radiated energy is shown as a function of the proper 
separation L between the apparent horizons. 
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For BL collisions, as explained in the previous section, there is a "favored" definition of bare mass. It is straight- 
forward in principle and easy in practice to compute a reasonable mass of each hole by going to large distances in 
the asymptotically flat region in which that hole is the only throat. For this definition of bare mass, Fig. [| shows the 
radiated energy as a function of L, with no anomalous behavior. This supports the argument that the low-L bump, 
for the Lindquist F — 0.75 curve in Fig. is due to an unphysical choice of bare mass. 



IV. RECOIL VELOCITY 



For unequal mass holes colliding along the z axis there will be momentum contained in the outgoing radiation due 
to the interaction of multipoles of different I. Moncrief [£3[ has considered a similar effect in the radiation emitted 
by a collapsing star. 

The general expression for the rate at which z momentum, is radiated is given by 



^1^J_V- / (*-!)(* + 3) d^ di/jt+t 



dt 16?r y (21 + l)(2l + 3) dt dt 

The mixture of I = 2 and 3 is dominant in the sum (^) when the holes start close but the next mixture becomes 
more and more important as they start farther and farther apart. As a result of radiation emission, the final hole 
formed will aquire a velocity (relative to the asymptotic frame in which the colliding holes were initially at rest). This 
recoil velocity is 

1 [ dP z 

v z = - — dt—. (14) 
M dt v ' 



In terms of (13) the recoil velocity is 




(15) 

Now the Zerilli functions ipz(t) and ipg+i(t) must be evolved numerically and the overlap integral in (|l^) must be 
numerically computed. The functions ipt(t) and ipe + i(t) have a fixed form and the details of the collision influence 
only the amplitude. The overlap integral in ( |l5|) therefore needs to be integrated only once. [See (|l2|).] 

This recoil velocity is always negative. Since our convention is always to put the more massive of the colliding holes 
on the positive z axis and the less massive on the negative side, this means that the final hole moves in the direction 
from which the smaller mass approached. Results are shown in Fig. ^| for the recoil velocity, in the case of BL initial 
data, as a function of separation. 



FIG. 6. Recoil velocity for BL holes. The recoil velocity of the final hole is shown as a function of initial separation for 
several ratios of bare mass. 



At the low values of L at which the CLAP should be reliable, the highest values of recoil velocity occur for F ~ 0.3. 
The recoil velocity results for ML initial conditions, shown in Fig. [7] are roughly similar. Mass ratios F ~ 0.1 — 0.3 
produce the highest recoil velocities. 

FIG. 7. For two ML holes the recoil velocity of the final hole formed is shown as a function of initial separation, for several 
different values of F, the ratio of Penrose bare masses. 

The figures both for ML and BL initial conditions include recoil velocities of many hundreds of km/sec, proper 
velocities that would be astrophysically interesting. The large values of v z correspond, of course, to large initial 
separations. Whether Figs. || and or [7] actually contain results of astrophysical interest depends on whether the CLAP 
fails at the separations which predict large recoils. This question then, gives us a very specific motivation for turning, 
in the next section to a consideration of the range of applicability of CLAP. 



G 



V. RANGE OF VALIDITY OF CLAP 



There is nothing inherent in the linear perturbation theory underlying the CLAP to indicate how small the expansion 
parameter e must be to have answers accurate to, say, 10%. In other words, we do not immediately know "how close 
is close enough." The only a priori justifiable way of finding the range of validity (aside of course from full numerical 
solutions) is carry out second-order perturbation theory, and to see at what e the first and second-order answers differ. 
The formalism for higher order computations has been developed fl(i| , and successfully applied to the equal mass case 
p| . Higher order perturbation theory, though much easier than full numerical solutions, is lengthy and tedious. For 
this reason it is useful to look for easy rough indices of validity of CLAP. 

To understand the importance of an index of validity one can look at the application of CLAP to the collision of 
equal mass ML holes [§]||. Here a comparison can be made with the results of full nonlinear numerical analysis. Both 
methods predict an increase of radiation with increasing initial separation. The computed energies agree reasonably 
well (within a factor of ~ 2) for initial separations out to around L/M ~ 4. The energy for infall from extremely 
large distances is larger only by a factor ~ 2 — 3. It may be that these features are generally true: the magnitudes 
of radiation quantities at the limit of CLAP validity may be within better than an order of magnitude agreement 
with the large separation limits. In that case we can make good estimates for the large separation case by taking the 
values at the CLAP limit. 



An easy index for this limit has been suggested by Suen |24|. One can take the nonlinear initial data, extract 
perturbation quantities and compute Moncrief's two even parity gauge- invariant functions q\ and q2- For linearized 
data the linearized Hamiltonian constraint gives 92 = 0, so the magnitude of q2 is an indication of "how nonlinear" the 
initial data are, and therefore, presumably, how accurate linearized computations are for the evolution. In Figs. ||, and 
^, we plot the value of 92/91 as a function of radius. To make this nonlinearity index more plausible we have modified 
it in two ways. The factor of (r/M — 2) corrects for the divergence of (72 at the horizon. If the mass ratio 7712/7712 
becomes very small, violations of sphericity must scale as 77i 2 . We would then find that violations of nonlinearity were 
very small regardless how far apart the holes start, and regardless how badly the CLAP fails. To take into account the 
scaling of nonlinearities as m 2 , we place a factor F in the denominator of the nonlinearity index. In the two figures 
we normalize separation with the ADM mass M, rather than with bare masses. Since the curves refer to a fixed set 
of masses, the ratio of M to mi + to 2 is fixed. 

The first figure shows the result for equal masses and ML initial geometry (see also j3j). In this case we know, 
from numerical relativity results, that the CLAP fails at around L/M ^3 — 4. The results in Fig. || show a dramatic 
increase in the hamiltonian violation as L/M increases beyond ~ 2 and this coincides, approximately, to the value of 
L/M for which CLAP starts to fail. 

In Fig. [9] we show the equivalent results for ML holes with a ratio of Penrose mass F — 0.1; in the previous section 
we saw that ratios roughly around this value maximize the recoil velocity. The results in this figure suggest that at 
L/M = 0.83, CLAP must fail badly. We know, from the particle limit study |^| that CLAP estimates are quite 
accurate (within 20%) up to around L/M ~ 1.5 This implies that the CLAP method can work well even when there 
is significant "linearity violation" in the initial data. It also implies, unfortunately, that there may be no easy reliable 
way to esitmate the validity of CLAP calculations. 

Let us now tentatively accept, from the particle limit results that for F = 0.1 the CLAP method is valid out to 
around L/M ~ 1, and that the energy and linear momentum at this limit are, to order of magnitude, as large as 
they would be for initially infinite separation (as in the case of equal mass holes). The plots in Figs. |^ and ^ then tell 
us that recoil velocities will always be well short of the several hundred km/sec values that would be astrophysically 
interesting. 

FIG. 8. The linearized constraint violation for two equal mass ML holes, F = 1 according with Penrose, is shown for several 
values of initial separation. Here r* is the coordinate r + 2M In (r/2M — 1) for the Schwarzschild geometry. 



FIG. 9. The linearized constraint violation for ML holes with a ratio F = 0.1 of the Penrose masses. 



VI. CONCLUSIONS 



The application of the CLAP method to collisions of unequal mass holes has provided one very tentative conclusion: 
The linear momentum generated in the collisions will be far too small to produce astrophysically interesting recoil 
velocity of the final hole. 



7 



The study has been much more successful in raising questions and in uncovering difficulties. The most important 
difficulty is the choice of the bare mass of the holes participating in a collision. It should be understood that this 
is not an issue specific to CLAP estimates. Rather, CLAP estimates were used to probe it. The general problem 
is: How do we best characterize an individual hole when it is interacting strongly with another gravitational source? 
And this issue is part of the broader question: How do we set up initial data to represent a black hole configuration? 

Astrophysical models, or Newtonian physics applied to the early dynamics, can give us an initial configuration in 
terms of simple physical parameters (mass of the holes, separation,...). To turn this into appropriate initial data for 
numerical relativity we must know the correct general relativistic interpretation of the classical picture. This will be 
of crucial importance to the use of numerical relativity to study three dimensional black hole collisions. Codes to 
evolve black holes tend to be unstable. It is therefore important to apply the codes "at the last possible moment," 
i.e., only to follow the last orbit or last few orbits of black hole coalescence. It will be necessary then for those codes 
to begin with initial data for holes which are already interacting fairly strongly and it is in just this case that there 
is the most ambiguity in translating the Newtonian concept of the mass of a hole, into a parameter of the relativistic 
initial data. 

We have considered three measures of bare mass. The "area-mass," and the "Lindquist mass" are well defined 
only for ML data, for which there are symmetric Einstein-Rosen bridges. The "Penrose" mass measure (though not 
generally applicable) can be computed for either ML or BL initial data (and in the latter case agrees with the obvious 
choice of bare mass in the asymptotically flat region of a throat). We have found very mixed messages about the 
physical meaningfulncss of the various mass measures when the holes are initially close. One somewhat surprising 
message is that the area-mass and the Lindquist mass of each throat are numerically quite close. They are even 
closer if we make a more subtle comparison: the ratio of area/Lindquist masses of two holes, and the ratio of Penrose 
masses. This agreement is not a manifestation of some general necessity of all mass measures to agree. The Penrose 
mass and the area/Lindquist mass are markedly different for close separations. 

The agreement of area-mass and Lindquist mass is unexpected, because they are based on such different criteria. 
The area-mass connects mass to the area of the minimum section of the Einstein-Rosen bridge, just as if it were an 
isolated Schwarzschild hole. It takes no direct account of the presence of a second throat. The Lindquist mass, by 
very sharp contrast, uses asymptotic masses of "images" used to form a single hole and subtracts an expression for 
binding energy due to gravitational interaction of the images. The binding energy is computed by Newtonian physics 
applied in the conformally flat space underlying the initial data. (Both the sum of asymptotic masses and the binding 
energy are divergent, but the sum is not.) It is often the case that when two very different ways of measuring a 
physical quantity agree it is taken as good evidence that the measurement is valid. This would suggest that we take 
seriously the area/Lindquist bare mass, perhaps for a wider class of problems. But the strange behavior seen in Fig. ^ 
for mass ratio F — 0.75 suggests even more strongly that this bare mass measure can be misleading. 



ACKNOWLEDGMENTS 



We would like to thank Carlos Lousto for many useful conversations, and also Madhavan Varadarajan for discussions 
of Penrose mass. We thank Carleton Detar for help with computational aspects of the work. ZA was supported by a 
PhD grant under PRAXIS XXI administered by JNICT (Portugal). This work was partially supported by NSF grant 
PHY-9507719. 



APPENDIX A: THE BRILL-LINDQUIST SOLUTION FOR TWO BLACK HOLES 

In this appendix we evaluate the coefficients in the expansion (JsJ) of the Brill-Lindquist solution (^). 
When \R\ > \Ri\ and \R\ > I-R2I the expansion of (||) in Legendre polynomials gives 

7 < = Mm + ( m>m ] - (A1) 

We can completely characterize the two holes by two parameters in the flat space: the "distance" zq = \R2 — Ri\ and 
the "mass ratio" C = a^/ai- Choosing the origin of coordinates in this space at the ficticious center of mass, we can 
rewrite 7^ as 

'-S^pTc?' 01 ^- 1 ^ (A2) 
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where we chose the convention that hole 1 is on the positive z axis, hole 2 on the negative side and < C < 1. Since 
the total mass of the spacetime at the moment of time symmetry is 

M = 2(a 1 + a 2 ), (A3) 



a i 



M 2(1 + C) 



(A4) 



we obtain 



^-^ l WTcW [cl + { - l)lc] - (A5) 

The two parameters zq and C do not have a direct physical meaning. We therefore introduce the ratio of the bare 
masses of each black hole, F = m 2 /m 1 , and the separation between the apparent horizons, L, measured along the 
axis of symmetry and we express z and C (and hence 7;) in terms of them. 

If we choose, for simplicity, the origin of coordinates to be at singularity l,thcn L, the proper distance between 
apparent horizons is 

Here z\ and z 2 are the z-axis intersections of the apparent horizons surrounding holes 1 and 2, respectively. To find 
zi , z 2 we numerically integrated the system of ODEs that determine all the extremal two dimensional surfaces of the 
BL solution (see Q for details), and we searched along the segment of the z axis between the positions of the two 
holes (i.e. between z = and z — z ) for the critical values z\ and z 2 at which the extremal surfaces are closed. The 
bare masses of hole 1 and hole 2 are jl3| : 

mi=2ai(l + — ) (A7) 

zo 

m 2 = 2a 2 (l + — ), (A8) 



and hence 



1 



M 



F=^= - 2(1+ c C) l] C. (A9) 



mi 



2(1+C) z 



Solving for C in terms of F and zq we obtain 



(F - 1)(2 + f ) + Ml Ff{2 + f f + 16F 

c = ¥ . (A10) 



Typically to characterize a BL solution, we choose an F and a set of zq values. Then from ( A10 ) we obtain the 
corresponding set of C values and from (Rq) the L value corresponding to each pair (C, zq). 



APPENDIX B: THE MISNER-LINDQUIST SOLUTION FOR TWO BLACK HOLES 

In this appendix we present the ML solution and some additional details relevant to the CLAP of this solution. We 
will follow closely Lindquist |L6) who derived such a solution in implicit form. Our notation is also the same as his 
with minor changes. 







1. The solution 



The conformal factor can be written in the explicit form (^) 



\X G r . 



(Bl) 



Where d\ is the position of sphere 1, of radius a and ei the position of sphere 2, of radius b, relative to a given 
origin in the flat space (the spheres in the flat space correspond to the throats in the three geometry) and d n are the 
positions of the image poles of d\ , e n are the positions of the image poles of e*i with a n and b n being the corresponding 
weights. These coefficients obey the recursion relations 

• If n is even 



d n = ei - 



ei + d n -i 



ei + dn-i 



= di- 



b n 



di + e n -i 



di + e Tl 



-b n - 



If n is odd (n > 3) 



d n = di- 
ai = a 



^ a 

bi 
b„ 



di + d n -i 



di + d n _i 
b 2 



0>n-l 



ei 



ei + e n _i 



= b 



d + e„_i 



whose solution is (see for example Smythe pq] ) 
• If n is even 



fea sinh(n + 2)/io + a 2 sinhn/io 

On = ei - c(l ) 

c z sinn n/io 

ab sinh 2/io 

a n = — r-j 

c sinn n[io 

6a sinh(n + 2)/io + fe 2 sinh n/io N 
e„ = rfi - c(l J—— ) 



c 2 sinh n/io 



&n = 



ab sinh 2/io 
c sinh n/zo 



(B2a) 
(B2b) 
(B2c) 
(B2d) 



• If n is odd 



di - c(l 



sinh(n + l)/io 



sinh(n + l)/i + f sinh(n — l)/io ' 



(B3a) 
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ab sinh 2/iq 



6sinh(n + 1)/xq + asinh(n — 1)/Xq 



(B3b) 



sinh(n + l)/Lto 
sinh(n + l)/io + sinh(n — l)/Lto ' 



e n = ei - c(l - -7-^ — — fc - — — — ) (B3c) 



b = ab sinh2/i d 

™ asinh(n + l)/io + &sinh(n — l)/io 

where c is the distance (in the flat space) between the centers of the two spheres 

c = d l +e 1 (B4) 

and /iQ is given by 

c 2 _ a 2 _ b 2 

cosh 2^o = 7. r ■ (B5) 

2ao 



2. The choice of parameters 

In the background three dimensional Euclidean space there are two natural dimcnsionless parameters: the ratio 
C = b/a of the radii a and b of the spheres 1 and 2, and the ratio of the distance between the two throats (i.e. the 
distance between the centers of the spheres c = d\ + ei) to the radius of one of the throats say a: D = c/a. We will 
always restrict attention to the cases < C < 1 i.e. hole 1 on the positive z axis is greater or equal to hole 2 on the 



negative z axis. We must also have D > 1 + C and from (B 1) 



1 D 2 — 1 — C 2 
Mo= 2 COsh_1 [ 2C ] (B6) 

These two parameters, like C and Zq/M in the BL solution, completely characterize the Misner-Lindquist solution. 
As in the BL solution, however, they do not have a direct physical meaning and we must look to other way of 
parameterizing the ML solution. One parameter which can be defined in a natural way in the ML solution, the 
distance between the two holes defined by the length of the geodesic which threads through their corresponding 
Einstein-Rosen bridges, L, is discussed in the first part of this section. Another physical parameter is the ratio 
between the "bare" masses of each ML hole, F = m%lm\. As discussed above, we lack a unique, well defined, notion 
of bare mass of each hole. The Lindquist masses presented on the second part are just a possibility. Once we have a 
prescription that gives the masses we can parametrize the ML solution either by F and L/M or by F and L/(mi+TO2)- 

Due to the difficulty of obtaining C and D as functions of these physical parameters, we will instead derive all the 
quantities in terms of C and D. 

As we will see the restriction D > 1 + C, equivalent to fj, > 0, is required in order to obtain convergence of the 
different series which appear in the next sections. 

a. Invariant distance of separation 

Instead of using the distance between the throats (i.e. the parameter D defined above) as measured in the Euclidean 
space, it makes more sense physically to use the invariant separation distance between the two holes defined as the 
length of the shortest closed path which threads through their corresponding E-R bridges. Such a curve is a geodesic 
of the initial slice. To evaluate its length it is convenient to work in bispherical coordinates instead of cartesian ones. 

In these coordinates the metric of the initial surface is 

ds 2 = $ 4 ( -J- ) 2 [dfi 2 + drj 2 + sin 2 i]dip 2 } . (B7) 

cosh [i — cos rj 

We introduce 

a = /csch/ii (B8a) 

di = f coth m (B8b) 

b = /csch fi2 (B8c) 

ei = /coth^ 2 (B8d) 
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where 



Mi + Ma = 2^o 

. i sinh 2/io 
/i2 = sinh 



/ = b sinh /i 2 . 



(B9) 
(BIO) 
(Bll) 



Substituting (|B8|) in (B2) and (B3) we obtain a convenient form for the coefficients of $ in (Bl). (Compare with (4.7) 
in Lindquist |l6|| when the two black holes are equal.) 



• If n is even 



a n = b n = /cschn/Jo 
d n = e n = / coth n^io 



(B12a) 
(B12b) 



• If n is odd 



a n = /csch [(n + 1)mo - Ma] 
6„ = /csch [(n - l)/x + Ma] 
rfn = /coth[(n + l)/x - Ma] 
e„ = / coth[(n - 1)mo + M2] 



(B13a) 
(B13b) 
(B13c) 
(B13d) 



and the conformal factor becomes 

$ = 



E 



1 



\/ cosh(/x + 4n/x ) - cos 77 
1 



y cosh[/i + 4n/io + 2/^2 ] — cos 77 



(cosh/i — cos 77) ■ 



(B14) 



which converges if and only if /Uq > 0. 

We note that this parametrization implies a specific choice of origin in the Euclidean space, however the geodesic 
distance between the two throats is the same independent of the choice of origin. It is thus just a matter of convenience. 
The equations of the throats x 2 + y 2 + (z — di) 2 = a 2 and x 2 + y 2 + (z + ei) 2 = b 2 become fj, = \i\ and /i = —^2 
respectively. 

The geodesic of interest is ip = 0, 77 = ir. Therefore the distance of separation between the holes is 



L = 



or explicitly 



/ / $W) 



M2 



1 



cosh /j + 1 



d[i 



(B15) 



l = 2/ a 



El 



4n/j 



sinh(/j2) sinh(2nMo) 
2n/Jo - M2 2n/j + M2 



sinh(2n/io — Ma) sinh(2n/io + M2) 
Again this series converges if and only if /iq > 0. 



(B16) 



b. The Lindquist mass 
The bare masses of the holes, according to Lindquist are 
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ii = 2^[a 2 „-i + 6 2 « + ^ 



a2n-l«2r? 



1 |<^2n-l — <^2m| 
Q2n-l02m-l , v%nfl>1m . &2n&2m-l 

\d2n-1 - e 2m _i| |e 2 „-d 2m | |e 2n -e 2m _i| 



(B17) 



"i 2 = 2 > la 2n + b 2n _i + > , — : + 



n=l 



02n&2 



m= l \d-2n — d2 m -l\ 
&2n-lG2m-l , &2n-l&2?i 



M2n-e 2m | |e 2 „_i - d 2m _i| |e 2 n-i e 2r 



which using the coefficients in ( B12 ) and ( |B13| ) become 

2 



mi 



= 2/$>{ 



•^-^ sinh2n^o sinh(2n/Xo — ^2) 



1 



sinh(2n/io + \li 



(B18) 



(B19) 



m 2 



2/E"i 



sinh[2(ra - l)/x + M2] 
Both series converge if and only if fiQ > 0. 



4 sinh2n/io sinh[2(?i + 1)/xq — /Z2] 



} 



(B20) 



3. The choice of origin 

In order to determine the energy and the recoil velocity we need to choose the origin in such a way that the dipole 
term in (@) vanishes. This requires that the origin be located at the "center of mass" (CM) of the "masses" a n , b n , 

00 

^ a n d nC M + b„e n cM = . (B21) 

71=1 

(The a n and b n do not depend on the choice of origin). One systematic way of determining the d n cM, e-nCM is the 
following 

• We choose an arbitrary origin, fo r exa mple t he o ne for which bispherical coordinates can be introduced and pick 
the corresponding d n and e„ in ( B12| ) and (B13). 



We next find the center of mass position using 



X)n=l ~ a 2nd2n + 0S2ra-l<fcn.-l + ^2nC 2 „ — 6 2 „_ie 2n _i 

ZCM = ^ — t • (B22) 

Z_-/n— 1 n ' n 

• Finally we determine the position of a given image relatively to the CM, using the formula 

zwm =Zi - z C M (B23) 

or in scalar form 

dinCM = -d 2n - z C m (B24a) 
e2n+iCM = -e 2 „+i - zcm (B24b) 
&2nCM = e 2 „ - zcm (B24c) 

d2n+lCM = ^2n+l — ZCM (B24d) 
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4. The 7^ coefficients 
Comparing (ji|) with (||), we obtain the coefficients 7^ 



v - zJif Her) + X? (ncr) ] ■ (B25) 



Af Af ' AT M 
n=l 



5. The mass of the system 



To determine the ji above we need to know a dimcnsionlcss parameter, say the radius of sphere 2, b, per unit of 
mass. This can easily be done by expressing the mass in terms of the coefficients in the conformal factor. The ADM 
mass of the spacetime at the moment of time symmetry is the coefficient of 1/r in an expansion of ds 2 in inverse 
powers of r. The result is 

M = 2j2 a n + K ■ (B26) 
n=l 



If the explicit expressions in (B2) and (B3) are used, this becomes 



00 9 1 

m = vY\- - 



J sinh 2n/xo sinh(2n/xo — /U2 
1 



} (B27) 



sinh[2(n - l)// + M2] 

from which we can express b in terms of C and D. We note that once more the series converges if and only if fiQ > 0. 

APPENDIX C: THE PENROSE QUASI-LOCAL MASS OF ONE ML BLACK HOLE 

In Penrose's approach a complex quantity A a p, the momentum-angular momentum twistor of the source inside a 
2-surface S (with topology S 2 ) is defined and the total mass threading through S is Jl7| 

m 2 P = -\A a0 A aP . (CI) 

Certain problems with this definition remain. In order to evaluate this mass we need a Hermitian "norm" defined 
for the surface S as a whole. This norm can only be defined unambiguosly if S can be embedded in some conformally 
flat space-time in such a way that both its intrinsic geometry (induced metric) and the quantities characterizing its 
extrinsic curvature are unaltered from their values when S is embedded in the true spacetime. However this is not 
a difficulty for the ML solution (and BL solution), since any time symmetric conformally fiat hypersurface can be 
embedded in a conformally flat 4-space as a surface of constant t with the same intrinsic and extrinsic geometries and 
hence also all 2-surfaces lying on this hypersurface. 

In fact Tod Jl8| has evaluated the Penrose quasi-local mass associated with a 2-surface S obeying the restriction 
above, including the time symmetric case. He showed that in the BL case, the Penrose mass of each hole coincides 
with the ADM mass and that for both BL and ML solutions the total Penrose mass of the spacetime at the instant 
of time symmetry is the same as the total ADM mass. He also derived the mass of each ML hole, when the masses 
are equal. 

In this appendix we use Tod's results to derive the Penrose mass of each ML black hole when their masses are 
unequal. We will use a signature (-) ). 

Let t denote the unit timelike vector normal to the time symmetric 3-surface E. The Penrose mass "enclosed" by 
the two surface S lying in E with normal n is 

ml = P a T - \X AB ^ab - ^ A ' B '~Pa'B' (C2) 
which we can rewrite in terms of spinors as ( see Tod p8|)) 
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in. 



I = -2t A , (B P*t c '( B P§ + P b t b -(PH b 



where 



-T AB ,A',C ,, , xA'B'.A .C'-n (ni\ 
A t B t A ,/J,AC + A t B ,t A HA'B' 

\iABt B A , = ^ / {2(rlV$)V$ - (V$) 2 ri}dS (C4) 
i Q P a = — / {-$n.V$ - 2(n.V$)(i?.V$) + 

27T 7 

(V$) 2 (H.R)}dS (C5) 
«A'(b-P# - ^ / (V$ x n) + (V$) 2 n xfi- 

2(V$.n)(V$ x (C6) 
2tt J 2 

ii? 2 (V$) 2 ) + i?(2(n.V$)(i?.V$) - (n.i?)(V$) 2 + 

$n.V$) - V^(rt.R) + R 2 (n.V^>)}dS (C7) 

with all expressions on the right hand side written in terms of a three vector notation on flat three-space. To evaluate 
the integrals we transform them to integrals over a volume spanned by S 

H AB t% = ~ J V$V 2 $d 3 a; (C8) 

t a P a = -— I ($ + 2R.V$)V 2 <S>d 3 x (C9) 
2vr J 

tA'( B P A) = \ J (R x V$)V 2 $d 3 2; (CIO) 

t% X A ' B ' = — [ [2R.V&R + $R - R 2 V$]V 2 <f>d 3 x . (Cll) 



2tt _ 

In the ML solution the conformal factor (|J) leads to 

v$ = y a „ dn ~ R +b n i n ~ R (ci2) 

and 

oo oo 

V 2 $ = -4tt[J2 a n S(R - d n ) + b n S(R - e„)] . (C13) 

n— 1 n— 1 

In order to evaluate the Penrose mass of the hole 1, we choose the 2-surface to be the sphere with center at R = d\ 
(the throat of radius a). Remembering that all the images located at R — d<i n —\ and R = em lie inside that sphere 
we obtain 



15 



oo oo 



, \ ^ \ ^ / d 2n — dn-\ , e 2n _i — d 2 i-i . 
i=i„=i |«2n - a2i-i| d |e 2 „-i - a2i-i| d 

oo oo T* 

4« > > 6 2j 6 2 „_i p =-,3 + a2«— ; — ) 

|e 2n -i-e 24 | 3 |rf 2 „-e 2i | 3 



oo oo oo , 2 _ ,2 

(rP a ) 1 = 2^ a 2,-i+6 24 + 2^^a 2l _ 1 [a2„^ rl 2 '-' 



= 1 n=l 



|d 2 i_i - d 2n \ 3 



p 2 - d 2 x x p 2 -p 2 d 2 p 2 

, e 2n-l a 2i-l i , p^^t r. c 2n-l e 2i , fl 2n e 2z ] 

\d 2i -i -e 2n -i\ 3 |e 2j -e 2 „-i| 3 \e 2l ~ d 2n \ 3 

[*A'(b-Pa)]i = 

since all images lie along the z axis, and 

00 00 00 »2 J* ^/2 J* 

(4A i,B ')i = -2i V om-!^-! + b 2i e 2i - 2i V V a 2 ,_ 1 [a 2 „ 2 " ~ ^-f 2 



i— 1 i— 1 n— 1 

2 J* j2 -* oo oo 2 -* 2 -* 
e 2n-l"2i-l — «2i-l e 2n-l i o ■ J, ri, e 2n-l e 2i — e 2i e 2n _i 
+ &2«-l ] - 2 * 2^ 6 2#2n-l p ~ ,3 

\d 2l -i -e 2 „_i| 3 |e 2j -e2„_i| 3 

+ a 2«— r ,3 J 

|e 2j - d 2 n| 3 

and the Penrose quasi-local mass of black hole 1 will be 

m\ P = l(P a t a )i] 2 + 2(\ AB ti) 1 {t c A , l jLAc)i ■ (C14) 
Similarly if we take for volume of integration the sphere of radius 6 we get 



oo oo oo ,2 _ j2 

(t«P ) 2 2^>, ■ fo 2l :-i 2££ a^^K- - " 1 



2—1 i—1 n— 1 



\d 2i -i - d 2n \ 3 



P 2 - d 2 oo oo 2 _ 2 .2 -2 

, , c 2n-l a 2i-l i V^Vl, rj, fc 2n-l c 2t , «2n _ e 2i 

+ 02n-l-^ - 2 > > »2i 02n-lT^; ~ TT + a 2' 



|rf 2 j_i - e 2n _i| 3 l= i„ = i |e 2i -e 2n _i| J |e 2l -d 2 n| 3 

[*A'(b-Pa)]2 = 



oo oo 



d\ n d 2l -\ - d^i-^n 



(tg,\ AB ) 2 = -2i^a 2 id 2 i + b 2 i-le 2 i-l + 2i'^2'^2a 2 i-l[a 2n - 
l=l i=ln=l M2i-l-rf2n| 3 

, e 2 „_ 1 c(2j-i - a2 i _ie2„-i 1 l0 .>r^\r^, r , e 2n _ 1 e 2 i - e 2t e 2n -i 

+o 2n -i + 2i > > b 2l [b 2n -i -; - r; 

\d 2i -i - e 2n -i\ 3 |e 24 -e 2 „-i| 3 

d\ r fi 2 i - e 2l d 2n 

+ a 2n—- J \ 

and the corresponding Penrose mass will be 

™ 2 2P = l(Pat a h} 2 + 2(\ AB 4') 2 (t%»Ach ■ (C15) 
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To compute the total mass at the instant of time symmetry we choose for surface of integration a 2-sphere surrounding 
both throats. Then each term above add up and we get 

(MAB^A')t = (MAB*A')i + (M-AB*A')2 = 

oo 

{t a P a ) T = {t a P a )i + (t a P a ) 2 =2j2ai + bi 

i=l 

\Pa'{bPa)]t = (*A'(B-PA))l + (*A'(B-Pa))2 = 

and 

OO 

(t%\ A ' B ') T = -2iJ2(aidi + hci) 



and hence 



M = 2^Tai + bi (C16) 



which coincides with the ADM mass (appendix B, expression ( B26[ )), As expected, this mass is not equal to the sum 



of the two individual masses. The difference is the binding energy 

E B = M - toip - m 2 p . (C17) 
In the special case in which the holes have equal mass, we obtain 

(t a P a )i = {t a P a ) 2 = \M (C18) 

(HABt A ,)i = -(nABt A ,) 2 = -il = - 4 *Zs2^ S i n h 2 (n + m) ^ 

oo . 

( ,A \A'B'\ 1 1 A \A'B'\ t o- 2 • i /V^ COSU n/i 

{t B ,\ )i = —(t B ,\ )2 — —iJ — —2ia sinn/io(> 2 

^ sinh n/xo 

OO OO . , 

^ y-* cosh n/io cosn rn^o . 
, , smh 2 (m + n) u 

n=lm=l > 



giving 



a result obtained by Tod fli 



m ip = m 2P = 3 m2 + 2JJ (C20) 



[1] A. A. Abramovici a/., Science 256, 325 (1992); K. S. Thorne, submitted to Proceedings of Snowmass 94 Summer Study 
on Particle and Nuclear Astrophysics and Cosmology, eds. W. W. Kolb and R. Peccei (World Scientific, Singapore). 

[2] Proceedings of the November 1994 meeting of the Grand Challenge Alliance to study black hole collisions may be obtained 
by contacting E. Seidel at NCSA (unpublished). 

[3] A. Abrahams and R. H. Price, Phys. Rev. D53, 1963 (1996). 

[4] R. H. Price and J. Pullin, Phys. Rev. Lett. 72, 3297 (1994). 



17 



P. Anninos, R. H. Price, J. Pullin, E. Scidcl and W.-M. Suen, Phys. Rev. D 52, 4462 (1995). 
A. Abrahams and R. H. Price, Phys. Rev. D53, 1972 (1996). 
A. M. Abrahams and G. B. Cook, Phys. Rev. D 50, R2364 (1994). 



J. Baker, A. Abrahams, P. Anninos, S. Brandt, R. H. Price, J. Pullin, E. Seidel, preprint |gr-qc/9608064| (1996) . 
H.-P. Nollert, private communication. 

R. Gleiser, C. O. Nicasio, R. H. Price, and J. Pullin, Class. Quantum Grav. 13 L117 (1996). 

R. Gleiser, C. O. Nicasio, R. H. Price, and J. Pullin, to appear in Phys. Rev. Lett. 

J. Bowen and J. W. York, Jr., Phys. Rev. D21, 2047 (1980). 

D. R. Brill and R. W. Lindquist, Phys. Rev. 131, 471 (1964). 

C. Misner, Ann. Phys. (N.Y.) 24, 102 (1963). 

C. Misner, Phys. Rev. 118, 1110 (1960). 

R. W. Lindquist, J. Math. Phys., 4, 938 (1963). 

R. Penrose, Proc. R. Soc. Lond. A 381, 53 (1982). 

K. P. Tod, Proc. R. Soc. Lond. A 388, 457 (1983). 

R. Penrose, in Global Riemannian Geometry, edited by T. J. Willmore and N. Hitchin (Ellis Horwood Limited, Chichester, 
England, 1984), p. 203. 

V. Moncrief, Ann. Phys. (NY) 88, 323 (1974). 
F. Zerilli, Phys. Rev. Lett. 24 737, (19 71). 



C. O. Lousto and R. H. Price, preprint ^r-qc/9609012; (1996). 
V. Moncrief, Ap. J. 238, 333 (1980). 
W.-M. Suen, private communication. 

William R. Smythe, Static and dynamic electricity (McGraw Hill, New York, 1968). 
N. T. Bishop, Gen. Rel. Grav. 16, 589 (1984). 



18 




0.2 ^ -L- -J 

0.0 1.0 2.0 

L/(m 1p +m 2P ) 



1.5 



1.0 



0.5 



0.0 



-0.5 



0.0 



2.0 



4.0 
L/M 



6.0 



8.0 



1.60 





0.0 1.0 2.0 3.0 4.0 5.0 




0.0 2.0 4.0 6.0 8.0 




0.0 2.0 4.0 6.0 8.0 



L/(m1+m2) 



\ 



\ 



L/M=1.5 

L/M=2.66 

L/M=2.88 

L/M=3.1 



\ 



-3.0 — 
—22.0 



-12.0 



-2.0 

r7M 



8.0 



-2.0 



a; 

c\j 
i 



-7.0 



\ 



\ 
\ 
\ 



\ 



L/M=0.43 
L/M=0.69 
L/M=0.83 



\ 



■12.0 
-22.0 



r7M 



-2.0 



